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Theories based on the coupling between spin fluctuations and fermionic 
quasiparticles are among the leading contenders to explain the origin of high- 
temperature superconductivity, but estimates of the strength of this interac- 
tion difl"er widely [Ij. Here we analyze the charge- and spin-excitation spec- 
tra determined by angle-resolved photoemission and inelastic neutron scatter- 
ing, respectively, on the same crystals of the high-temperature superconductor 
YBa2Cu306.6- We show that a self-consistent description of both spectra can be 
obtained by adjusting a single parameter, the spin-fermion coupling constant. 
In particular, we flnd a quantitative link between two spectral features that have 
been established as universal for the cuprates, namely high-energy spin excita- 
tions in HI SI El El [7] and "kinks" in the fermionic band dispersions along the 
nodal direction [HI [9l llOj, IIH, 112] . The superconducting transition temperature 
computed with this coupling constant exceeds 150 K, demonstrating that spin 
fluctuations have suflicient strength to mediate high-temperature superconduc- 
tivity. 

Looking back at conventional superconductors, the most convincing demonstration of 
the electron-phonon interaction as the source of electron pairing was based on the quan- 
titative correspondence between features in the electronic tunneling conductance and the 
phonon spectrum measured by inelastic neutron scattering (INS) [THl [H]. The rigorous 
comparison of fermionic and bosonic spectra was made possible by the Eliashberg theory 
[15] which allowed one to derive the tunneling conductance from the experimentally de- 
termined phonon spectrum. Various difficulties have impeded a similar approach to the 
origin of high-temperature superconductivity. First, the (i-wave pairing state found in these 
materials implies a strongly momentum-dependent pairing interaction. A more elaborate 
analysis based on data from momentum-resolved experimental techniques such as INS and 
angle-resolved photoemission spectroscopy (ARPES) is thus required. These methods, in 
turn, impose conflicting constraints on the materials. In order to avoid surface-related 
problems, most ARPES experiments have been performed on the electrically neutral BiO 
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cleavage plane in Bi2Sr2Ca„_iCu„02(n+2)+5 [!]• However, as a consequence of electronic 
inhomogeneity this family of materials exhibits broad INS spectra that greatly complicate 
a quantitative comparison with ARPES data [7J. Conversely, compounds with sharp spin 
excitations, including YBa2Cu306+x, have generated problematic ARPES spectra due to 
polar surfaces with charge distributions different from the bulk [8]. Finally, an analytically 
rigorous treatment of the spin-fluctuation- mediated pairing interaction is difficult, because 
small expansion parameters used in the traditional Eliashberg theory (such as the ratio of 
Debye and Fermi energies) are missing [TB]. Because of these difficulties, widely different 
values have been quoted for the spin-fermion coupling constant [1]. 

The analysis of YBa2Cu306.6 data reported here was made possible by recent advances on 
several fronts. First, INS experiments on this material now consistently yield high-quality 
spin excitation spectra over a wide energy and momentum range [21 H]. Second, recent 
ARPES experiments on YBa2Cu306+a; [HI [12] have overcome problems related to polar 
surfaces and allowed the observation of superconducting gaps and band renormalization 
effects ( "kinks" ) akin to those previously reported in La- and Bi-based cuprates [8J . Third, 
calculations based on the two-dimensional (2D) Hubbard model have demonstrated Fermi 
surfaces, single-particle spectral weights, antiferromagnetic spin correlations, and d^^-yi 
pairing correlations in qualitative agreement with experimental measurements p!7l [T8l [T9l 
[20]. Numerically accurate solutions of this model can thus serve as a valuable guideline for 
a treatment of the spin fluctuation interaction in the cuprates. This is the approach we take 
here. 

Recent quantum Monte Carlo calculations of the 2D Hubbard model within the dynamical 
cluster approximation (DCA-QMC) |20] have shown that the effective pairing interaction 
can be parameterized in terms of the numerically computed spin susceptibility xiQ^ ^) 
the form 

Veff{Q,n) = ^u'xiQ,^), (1) 

where f7 is the coupling strength, and that this interaction generates reasonable values 
for the superconducting transition temperature Tc. Here we follow a similar strategy, but 
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use determined by INS, on high-quality detwinned YBa2Cu306.6 single crystals 

described previously In order to serve as input for the numerical calculations, we have 
used an analytic form of Imx that provides an excellent description of the INS data (see 
the Supplementary Information). Fig. 1 shows a plot of this form in absolute units. In the 
superconducting state, the spin excitations exhibit the well-known "hour glass" dispersion, 
with a neck at the wave vector Q = (tt, tt) characteristic of antiferromagnetism in the copper 
oxide planes and the "resonance" energy Q = 38.5 meV. (We use a notation in which the 
lattice parameter a and the Planck constant h are set to unity.) The lower branch of the 
hour glass appears to be influenced by materials-specific details. For instance, recent INS 
work on La2_2.Si'xCu04 indicates two characteristic energies [3, E], rather than the single 
resonance found in YBa2Cu306+x- The upper branch of high-energy spin excitations, on 
the other hand, is common to all copper oxides thus far investigated by INS in this energy 
range [21 |3l HI O EJ [TJ |2T]. Moreover, while the resonance in YBa2Cu306.6 disappears above 
Tc [1], the intensity of the spin excitations above ~ 50 meV is not noticeably affected by 
the superconducting transition and only decreases slowly upon further heating [21 S]. 

We extract the second parameter in Eq. 1, the coupling strength U, from a combined 
analysis of the INS data parameterized in this way and the fermionic band dispersions 
observed by ARPES on the same crystals (see Methods). As noted before, bonding and 
antibonding combinations of electronic states on the two Cu-0 layers in the YBa2Cu306.6 
unit cell give rise to two distinct Fermi surfaces (Fig. 2). The most prominent signature of 
many-body effects in the ARPES data, namely the "kink" along the nodal direction (cut 
#1 in Figs. 3), is highlighted in Fig. 4, where the bonding band is singled out by a proper 
choice of excitation energy. 

We now proceed to a quantitative analysis of the renormalization of the nodal band 
dispersion by spin fluctuations. Before describing the results, we take a look at the kine- 
matics of spin fluctuation scattering near the nodal points, where complications from the 
superconducting gap are absent. The spin fluctuations shown in Fig. 1 scatter electrons be- 
tween bonding and antibonding bands, as indicated by factors in the INS [21 H] and ARPES 
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[TT| [T2] cross sections. (Weak high-energy excitations corresponding to intraband scattering 
[22] are neglected here.) An analysis of our numerical results below shows that the scattering 
probability for electrons near the nodal points is greatly enhanced when energy-momentum 
conservation allows interband scattering into opposite nodal regions (green arrow in Fig. 2). 
A look at the INS data (green line in Fig. 1) reveals that this condition is satisfied by spin 
fluctuations of energy ~ 80 meV on the upper, universal, weakly temperature dependent 
branch of the hour glass. At this characteristic energy in the temperature range studied 
here, we therefore expect a weakly temperature dependent anomaly in the band dispersion, 
as experimentally observed. 

A self-consistent numerical procedure with a single adjustable parameter, the coupling 
strength U, was developed to quantitatively assess the influence of the spin-fluctuation in- 
teraction on the spectral function determined by ARPES (see the Methods section). Fig. 4 
shows that an excellent description of the nodal band dispersion over a wide energy range 
is obtained with U = 1.59 eV, in rough agreement with values found in earlier calculations 
based on phenomenological models of the spin susceptibility [Ij. In particular, both theo- 
retical and experimental results show deviations from linear behavior ("kink") for u > 80 
meV (arrow in Fig. 4). The corresponding mass renormalization at the nodal point is 
Re Za = 3.7. 

Fig. 3 shows a comparison of the calculated spectral weight to the ARPES intensity for 
all three cuts in Fig. 2. It is evident that the calculation yields an excellent description of 
the ARPES data set over the entire Brillouin zone without further fitting parameters. In 
particular, the low-intensity region ("dip") below the renormalized band in cut #3 can be 
understood as a consequence of coupling to the magnetic resonance at the neck of the hour- 
glass dispersion. As noted before [H [231 [21], the resonance wave vector (red lines in Fig. 1 
and 2) connects antinodal regions in bonding and antibonding bands, and the resonance and 
gap energies add up to the dip energy ~ 65 meV. The only noticeable difference between 
the numerical and experimental data is the width of the momentum distribution curves, 
which is substantially larger in the ARPES data, presumably at least in part due to residual 
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surface inhomogeneities [25]. 

Encouraged by the self-consistent description of INS and ARPES data, we proceed to a 
calculation of the critical temperature of the d-wave superconducting state arising from the 
exchange of spin fluctuations. A recent DCA-QMC study has shown that a good estimate 
of Tc can be obtained by using the same effective interaction as in the calculation of the 
self-energy [20j. For the set of parameters found above, the linearized gap equations (see 
Methods) yield the d-wave eigenvalue = 1.39 in the normal state (T = 70 K), correspond- 
ing to a transition temperature Tc = 174 K. In principle, the INS and ARPES spectra would 
now have to be remeasured at this higher temperature, the calculation repeated, etc., until 
self-consistency is achieved. However, as the spectral weight rearrangement of spin excita- 
tions in this temperature range is largely confined to low excitation energies, our estimates 
for U and Tc are not expected to change substantially (see Methods). In this context it 
is instructive to compare the eigenvalue at T = 70 K with the one obtained from the INS 
spectrum at 5 K that includes the "resonance", = 1.49. The enhanced eigenvalue implies 
that the redistribution of spectral weight of the spin excitations below Tc leads to an increase 
of the effective pairing strength. This lends support to an interpretation of the magnetic 
resonance and associated antinodal dip in terms of a superconductivity-induced feedback 
effect on the spin fluctuation spectrum [Ij. It is also consistent with the large 2Ao/ ksTc 
ratio. 

In summary, we have shown that data from two momentum-resolved experimental probes 
of a cuprate superconductor can be related in a quantitative fashion, in close analogy to 
the traditional analysis of the electron-phonon interaction in conventional superconductors. 
Our analysis is in overall agreement with conclusions drawn from prior work based on phe- 
nomeno logical spin excitation spectra and/or data from probes without momentum resolu- 
tion [H [23|, [2U |26l [271 12H] , and it resolves some problems that appeared in the context of 
these studies. In particular, models that attribute the nodal kink in ARPES either directly 
to the magnetic resonance [231 122] or to incoherent scattering processes from a node into 
gapped states at the antinode [21] generally predict that the kink is strongly modified by the 
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onset of superconductivity, whereas the experiments indicate at most a weak effect at Tc [lO] ■ 
The node-to-node interband scattering process mediated by weakly temperature dependent, 
universal, incommensurate spin excitations we have identified provides a straightforward 
explanation of this observation. As the incommensurate spin excitations persist into the 
optimally doped [Sj [7] and overdoped [6] regimes, it also explains the persistence of both the 
kink and superconductivity at high doping, where feedback effects related to the magnetic 
resonance are progressively reduced |27|. There is thus no need to invoke phonon scattering 
at this level [S]. While some contribution of phonons to the nodal kink cannot be ruled out, 
recent work has shown that it is hard to obtain a quantitative description of the kink based 
on the electron-phonon interaction alone [2H EQl El] . 

It was previously shown [22] that the change in the magnetic exchange energy between 
the normal and superconducting states is more than enough to account for the cuprate su- 
perconducting condensation energy. However, the crucial question whether the exchange of 
magnetic spin-fluctuations actually has the strength to give rise to high-T^ pairing, was not 
answered. Here we have shown that this interaction can generate d-wave superconducting 
states with transition temperatures comparable to the maximum observed in the cuprates. 
In any given material, especially underdoped cuprates such as YBa2Cu306.6, a variety of ef- 
fects not considered in our analysis can reduce the actual Tc, including vertex corrections of 
the spin fluctuation interaction (TU] , phase fluctuations of the order parameter, competition 
with other types of order, and pair breaking by phonons and impurities. It is also possible 
that phonons P, [30l El] or higher-energy excitations [33j contribute to the pairing interac- 
tion. However, our analysis indicates that the exchange of spin excitations already directly 
observed by INS is a major factor driving the high temperature superconducting state in 
the cuprates. 
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Methods 



The ARPES measurements were performed on the same YBa2Cu306.6 crystals used for 
the INS experiment, thus avoiding systematic uncertainties invariably associated with mea- 
surements on different materials. The details of the ARPES experiments have been described 
elsewhere [HI [I2] . Usually, YBCO single crystals cleave between the Ba-0 and Cu-0 chain 
layers, resulting in an effective overdoping of the Cu-0 layer closest to the surface [HI [T21 El] . 
A recent comprehensive study has revealed, however, that in some cases the ARPES spectra 
are dominated by a signal from the nominally doped Cu-0 plane [12] • Here we present data 
taken on a particular spot on the surface after one such successful cleave (Figs. 2, 3). The 
strong many-body renormalization of the band structure typical for underdoped cuprates 
(Fig. 4) as well as the anisotropic superconducting gap (cuts #2, 3 in Fig. 3) in the ARPES 
spectra demonstrate that contributions from the overdoped surface-related component are 
negligible. The superconducting component we observe corresponds to the nominal doping, 
as follows from the size of the gap and the temperature evolution of the coherence peaks 
which disappear above the bulk SC transition temperature [T2] . 

The self-consistent numerical calculation we have used is based on the self-energy diagram 
displayed in Fig. 4. The Green's functions G{k,uj) on antibonding {A) and bonding (B) 
bands can be written as [IHj 

GABik, u) = ^'^ + ''^'"^ (2) 

{ujZA,By - {it'"") -{ReZA,B{k,uJ = 0)Aky 

where is the superconducting gap, which we assume to be of the d-wave form = 
Ao(cosA;a; — cosky)/2 with Aq = 30 meV, and e^'"^ = e^'^ -|- ^a,b is the renormalized band 
structure. The unrenormalized band dispersions e^.' were derived from tight-binding fits 
to the ARPES Fermi surface in combination with additional information from band struc- 
ture calculations (see the Supplementary Information). We have found that the results 
of our calculations are quite robust against modifications of ef'^ (see the Supplementary 
Information). Finally, ujZA,B{k,Oij) = uj — ^ (YjA,B{k,uj) — S^^(fc, —uj)^ + iT^i is the mass 
renormalization function and ^A,B{k,uj) = ^ (TjA,B{k,uj) + ^A^^ik^ —^)^ the energy shift 
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function. Apart from an elastic scattering rate Fei ~ 30 meV, which accounts for impurity 
scattering, the mass renormahzation function is determined by the imaginary part of the 
electron self-energy I^a,b, which can be written as 

ImEA,B(A;,w) = — 5] / dn[n{n) + f{n-uj)]lmVeff{Q,n)ImGBAk-Q:i^-^) ■ (3) 

TTiV Q J-oo 

Here J2q denotes a sum of the in-plane momenta over the full Brillouin zone, n and / 
are the Bose- and Fermi-functions, respectively, and is the spin fluctuation interaction 
(Eq. 1). The real parts of T,a,b that enter into Eq. 2 are obtained by Kramers-Kronig 
transformations. Note that the self-energy in the antibonding band is determined by the 
interaction with the bonding band and vice versa. 

Together with Eq. 1, this defines a self-consistent system of equations with a single 
adjustable parameter, the coupling strength U. Starting with noninteracting values for the 
Green's functions, these equations were solved iteratively until convergence is achieved. The 
renormalized band dispersion and spectral weight, f{uj)lmG{k,uj), can then be compared 
to ARPES data. 

The linearized gap equations 

Xa Im <l>A,B{k,uj) = ^ E / dw'inicu - cu') + f{-u;')] (4) 



Im V^//(^ — k'^uj— uj') Im 



were solved for the same set of parameters. Note that the INS data used for the calculations 
were taken at T = 5 and 70 K, while the ARPES data were taken at 30 K. Since the changes 
in the superconducting gap and INS spectrum between 5 and 30 K are negligible, we use 
the 5 K INS results along with the 30 K ARPES data to determine the coupling constant 
U. As T is raised further, the superconducting gap decreases, and there is a shift of Imx to 
lower frequencies. However, we expect that U is unchanged for the range of temperatures 
of interest, because it is determined by weakly T-dependent high-energy processes. 

For discussions of the influence of a high-energy cutoff in x(Q, ^) and of a normal-state 
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pseudogap on the results of the numerical calculations, see the Supplementary Information. 
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FIG. 1: Intensity of spin excitations along Q = q{7r,TT) resulting from numerical fits to the INS 
spectra of YBa2Cu306.6 at T = 5 K. The maximum value corresponds to 600 /i^/eV/f.u. The 
green (red) lines mark wave vectors connecting nodal (antinodal) regions on different Fermi surfaces 
(Fig. 2). 
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FIG. 2: Results of tight-binding fits of the Fermi surfaces determined by ARPES for the bonding 
(soUd hne) and antibonding (dotted hne) bands. The blue lines denote the three cuts along which 
experimental and theoretical spectral functions are compared (Fig. 3). The green and red arrows 
indicates spin-fluctuation-mediated scattering processes discussed in the text. 
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FIG. 3: a-c ARPES intensities at 30 K along the three cuts shown in Fig. 2, obtained with photon 
energies 60 eV (cut #1) and 55 eV (cuts #2, 3). d-f Theoretical spectral functions at T = 5 K, 
multiplied by the Fermi function at 30 K (sec Methods). Intensities in bonding and antibonding 
bands in panel d were multiplied by 0.75 and 0.25, respectively, to take account of matrix element 
effects. The black lines in panels e and f show the unrenormalized antibonding band dispersion. 
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FIG. 4: Nodal dispersion of the ARPES data (crosses) for the bonding band compared with the 
same nodal dispersion of the calculation (solid line). Here, the coupling constant U was adjusted 
such that the slopes of the curves (renormalized Fermi velocity) match near zero energy. The inset 

shows a sclf-cncrgy diagram for the exchange of spin fluctuations. The wiggly line represents the 
interaction given by Eq. 1, and the straight double line is the self-consistent single-particle Green's 
function. 
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Fit to INS data 



The INS results were fitted by simple analytic forms which could then be numerically 
integrated in the self-energy and Tc calculations. In constructing these fits, the energy range 
was divided into upper {u) and lower {£) branches corresponding to energies above and below 
the resonance energy cUr- 

Im x(Q, ^) = Im Xe(Q, ^) + Im Xu(Q, (5) 

The data in the superconducting state at T = 5 K were fitted as follows. On the lower 
branch 

Here the frequency Q is measured in meV with the resonance frequency Qr = 38.5 meV. 
The prefactors are given by 

^'(") = l- l + exp[(L24)/3| « 

and 

= l+oxp|n-42] 

These factors essentially cut off the lower branch at frequencies above 42 meV and below 24 
meV (spin gap) . The momentum dispersion of the resonance mode is given by 

= + SaiQa - Qoaf + Sb{Qb - Qobf (9) 

with Qo = (<5oa) Qob) the antiferromagnetic wave vector, = 280 meV and — 455 
meV A^. The damping width 7q has an angular dependence of the form 

. A {Qa Qoa) /1A\ 

7q = 7a + ^7^,^, — (10) 
|v4 ~ Mo I 

where the index a denotes the a-axis direction in the Brillouin zone and b the 6-axis direction, 
respectively. The parameter values are ja = 1-5 meV and A7 = 6.5 meV. 
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For the upper branch we have 

Im XuiQ, ^) = NoNoMm<iipimrotiQ) _ ^ ^^^p^, (11) 
Here, F = 11 meV and 

= + [sl'^iQa - Qoaf + Sl'\Q, - Q^,fY'^ (12) 

where p 4 is the power of the dispersion of the upper branch (somewhat steeper than 
quadratic) with Sa = 4830 meV and Si, = 10065 meV A*'. The prefactor Nq^ is 1.35 and 

NJ9) = 1 — - (13) 

' l + expp-36)/1.5] ^ ' 

provides a cut-off of the upper branch for energies below 36 meV. The prefactor 

Naij>{n) = 1-0.2 exp [- (Q - Q^^^f / (24^)] (14) 

with ^l^ip — 47 meV and (Tdip — 2 meV accounts for the reduction in intensity just above the 
resonance. The prefactor N^ot has a momentum dependence of the form 

NrotiQ) = 1-0.3 |^ (Q"-^0^ J-^^^^-Wy (15) 

accounting for the 45 degree rotation of the signal above the resonance [SI, S2]. The overall 
normalization factor A^o is chosen such that the momentum integrated Im x{Qi ^) has a peak 
value of 16 n^/eV /i.u. in the superconducting state at 5 K, which is known from previous 
work [SI, S3, S4]. In order to account for lattice symmetry, the function Im x(^, Q) is cut 
at the zone center and continued periodically. 

Let us now discuss the fitting formula for the normal state data at 70 K. Again, the 
fitting formula is subdivided into an upper and lower branch. For the lower branch we have 



Imx^(Q,^^) = NoNoeNiinmm{n)NLm (16) 
{exp[-41n2[((g„ - Qia)/aaf + {{Qb - Qib)/(7bf]] 

+ exp[-41n2[((Q„ - Q2a)/<Taf + {{Qb - Q2b)/<Tbf]]} 

3 



Here, the two Gaussians describe an energy-independent incommensurability with Ua — 
0.125 r.l.u., (Tfe = 0.17 r.l.u., Qi = {Qxa^Qxh) = (0.575,0.5) r.l.u., and Q2 = {Q2a,Q2b) = 
(0.425,0.5) r.l.u. The normalization factors are given by A^o^ = 0.001807, 

^'(")-l- l + exp[(n-29.5)/2] ^"^ 

and 

^'("' = l+exp[(n-37)/2] - 

These factors essentially cut off the lower branch at frequencies above 37 meV and partially 

below 29.5 meV. The factor 

1 for Q > 27meV 



Nun m 



(19) 

n/27 for n < 27meV 



represents a simple linear decrease below 27 meV consistent with previous work [S5]. For 
the upper branch we have the same expression as in the superconducting state at 5 K except 
that 

Ndipifl) = 1 (20) 

and 

1 02 

^''") = '- l+exp[(0-37)/2| - 

Figure SI shows the momentum integrated Im x{q, ^) in absolute units obtained from 
the two formulae. The "sum-rule" integral 

for the two fits gives S = 0.071/f.u. for the normal phase (T = 70 K) and 0.070/f.u. for the 
superconducting phase (T = 5 K). These numbers are in reasonable agreement with each 
other. 

Figure S2 shows typical measured scans compared with our fit formula convoluted with 
the instrumental resolution function and demonstrates the good agreement. 

In order to check the sensitivity of our results to the high energy part of the magnetic 
excitation spectrum, we have repeated our calculations with a spectrum cut off at 200 meV. 



As a result we found that the couphng constant U had to be increased by 3 percent. Tc 
was found to decrease by 1 Kelvin, showing that the details of the high energy part of the 
spectrum do not affect very much. 

Further on, we have phenomenologically added a normal state d-wave pseudogap of 30 
meV into our calculation and found that Tc increases by 20 percent, due to the suppression 
of pairbreaking low-energy magnetic excitations. We conclude that the influence of the 
pseudogap does not alter our two major findings: the high value of Tc and the nodal kink 
generated by the upper branch of the hourglass. 

Fit to fermionic band dispersions 

As a starting point for our theoretical calculation we also need the unrenormalized dis- 
persions e^'^ for the bonding (B) and antibonding (A) bands of the two-layer system. In 
contrast to previous calculations, here we keep the renormalizcd Fermi surface fixed during 
the iterative solution of Eqs. (1-3) of the main manuscript, because the Fermi surface is 
known from the ARPES data. In order to achieve this, we keep the renormalized quantity 
e^'^ = e^'^ + Re ^A,B{k, to = 0) fixed during the calculation and obtain it from tight-binding 
fits to the ARPES data of the form 

I^A,B — 2tA,B(cos kx + COS ky) + 4t^^^ COS kx cos ky — 2t'^_^(cos 2kx + cos 2ky) (23) 

with the parameters 

HA = 556 meV tA = 409 meV t'^ = 150 meV t'X = 40 meV 
Hb = 417 meV te = 550 meV t'^ = 231 meV t'^ = 67 meV 

which yield excellent descriptions of the ARPES Fermi surfaces (Fig. 2). These pa- 
rameters are scaled such that the unrenormalized Fermi velocity at the nodal point in the 
antibonding band equals 5 eVA, as found in ab-initio bandstructure calculations [S6]. 

In order to test the sensitivity of the numerical results to the assumed unrenormalized 
Fermi velocity, we have repeated the calculations with an unrenormalized Fermi velocity of 



5 



4 eVA, a value that deviates from the ab-initio predictions for YBa2Cu306+a; but is close 
to the one found in Bi2Sr2CaiCu208+<5 [S7]. In order to reproduce the renormalized Fermi 
velocity, the coupling constant U has to be reduced from 1.59 to 1.23 eV, and the nodal 
mass renormahzation drops from 3.7 to 3.0. As shown in Fig. S3, the agreement of the 
numerical and ARPES data for the nodal dispersion at higher energies is worse than that 
obtained with the more realistic parameters in the main text. Nonetheless, the estimate of 
the critical temperature for d-wave superconductivity remains high. For the INS spectrum 
at 70 K, we obtain a = 1.28, corresponding to = 140 K. 

Correspondence and requests for materials should be addressed to B.K. 
(b.keimer@fkf.mpg.de) . 
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FIG. 1: Momentum-integrated spin excitation spectra at T = 5 K and T = 70 K according to the 
fitting formulae. 
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FIG. 2: Comparison of typical constant-energy scans at T = 5 K (black squares) with our fit 
formula convoluted with the instrumental resolution function (open cirles). The black lines are 
guides to the eye for the measured data points. High- and low-energy scans were measured at 
different kf, thus the intensities are not comparable. 
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FIG. 3: Comparison of the nodal dispersion measured by ARPES (crosses) and evaluated theoret- 
ically for unrenormalized Fermi velocities of 4 eVA (solid line) and 5 eVA (dashed line). 
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